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ri Abstract 

Using the AdS /CFT correspondence, we calculate the transport coefficients of a strongly 
^ interacting system with a non-abelian SU(2) global symmetry near a second order phase 

transition. From the behavior of the poles in the Green's functions near the phase tran- 
sition, we determine analytically the speed of second sound, the conductivity, and diffu- 
sion constants. We discuss similarities and differences between this and other systems 
with vector order parameters such as p-wacve superconductors and liquid helium-3. 
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1 Introduction 



A new and interesting holographic perspective on the physics of superfluids and supercon- 
ductors was provided by refs. [U |2]. These papers, which rely on the AdS/CFT correspon- 
dence [31 m |5] , provide a dual description of the superconducting phase transition as the 
instability of a charged black hole to develop scalar hair. Recalling that the AdS/CFT 
correspondence maps a strongly interacting field theory to a classical gravity description, 
this new perspective holds promise for deepening our understanding of superconductivity 
in strongly interacting regimes where BCS theory |6j is inadequatej^ 

This paper takes as its starting point the holographic model proposed in refs. [H [9| 110] . 
While in refs. [U [2] the black hole develops scalar hair at the phase transition, in refs. 
[H ini [To] the black hole develops non-abelian hair. More precisely, refs. [H [2] begin with 
gravity plus an abelian gauge field and charged scalar, while refs. [H |9l [10] omit the scalar 
and promote the abelian field to a non-abelian SU(2) gauge field. Recall that the AdS/CFT 
dictionary maps gauge fields on the gravity side of the duality to global symmetries in the 
field theory. 

We are intentionally vague about the distinction between a superfluid and supercon- 
ductor. As emphasized in ref. |llj . technically, the field theory dual to the black hole 
construction undergoes a superfluid phase transition, i.e. spontaneous symmetry breaking 
of a global symmetry. To interpret this transition as a superconducting phase transition, 
the global symmetry must be weakly gauged. For many questions about superconductivity, 
the distinction is irrelevant, and from the current-current two point functions we calculate 
below, we can extract meaningful conductivities. 

This SU(2) system has many desirable features compared to the scalar system. The 
scalar system appears to be less universal; one must specify a potential for the scalar, or 
at the very least a mass term. In comparison, up to the strength of the gravitational and 
gauge couplings, the form of the Lagrangian for the SU(2) system is completely specified 
by gauge invariance. Also, it appears more straightforward to embed the SU(2) system in 
string theory or other UV completions. (See however ref. |[12J for recent progress with the 
scalar system.) This embedding can move the SU(2) system out of the toy model realm 
and give us a microscopic Lagrangian for the field theory. Refs. |13| [H] discuss one possible 
embedding where the SU(2) is a global flavor symmetry group for = 4 super Yang-Mills 
broken to AA = 2 by the addition of two hypermultipletsj^ 

Our reason for choosing this SU(2) system is very simple: We can get analytic results 
near the phase transition. The differential equations that describe these holographic systems 
are nonlinear, and analytic solutions do not appear to be available in most cases. The 
studies mentioned above make extensive use of numerics to see the phase transition, to 



^See ref. [3 for a review of the limits of BCS theory when confronted with high temperature supercon- 
ductivity. 

^The most naive embedding does not actuaUy work. The original AdS/CFT correspondence gives us 
a duality between A/" = 4 super Yang Mills and string theory in the curved background AdS^ x . This 
string theory can be approximated at low energies by a gauged supergravity in AdS;,. The SU(4) gauge field, 
which maps to the SU(4) global R-symmetry in the field theory, has an SU(2) subgroup. Unfortunately, 
supersymmetry constrains the relationship between the gauge field coupling and the gravitational coupling 
in this model, and the gauge field coupling is too weak for a superconducting phase transition to occur [TD] . 
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calculate the conductivities and critical exponents. Analytic results, for example the low 
temperature approximation of the conductivity in ref. |2j , are scarce Our starting point 
is the remarkable observation in ref. [13] that the zero mode for the phase transition for an 
SU(2) gauge field in AdS^ has a simple analytic form. From this zero mode, we are able to 
extract a long list of properties near the phase transition: 

1. The speed of second sound near the phase transition. 

2. That the phase transition is second order. 

3. The conductivity and in particular the residue of the pole in the imaginary part of 
the conductivity. 

4. The system satisfies a London type equation that implies a Meissner effect. 

5. A large selection of current-current Green's functions in the hydrodynamic limit, and 
that they satisfy the appropriate non-abelian Ward identities. 

To clarify the title of the paper, recall that in a two component fluid, there are typically 
two propagating collective modes. The first mode corresponds to ordinary sound in which 
the two components move in phase. The second mode corresponds to second sound in 
which the two components move out of phase. Typically, ordinary sound can be produced 
by pressure oscillations while second sound couples much more strongly to temperature 
oscillations [17] . 

The order parameter for the phase transition in our SU(2) model is the set of non- 
abelian global SU(2) currents, ja. As pioneered in ref. [8], we introduce by hand a chemical 
potential in the third isospin direction which induces a charge density, (jg) ^ 0, that breaks 
both the global SU(2) symmetry to a U(l) sugroup and also Lorentz invariance. There 
is a superconducting phase transition at a critical temperature Tc, below which a current 
develops orthogonal to the third isospin direction that completely breaks the residual U(l) 
symmetry and also breaks the remaining rotational symmetry of the system to U(l). For 
convenience, we take this current to be in the direction (j'^J), leaving a rotational symmetry 
in the yz-plane. 

The fact that rotational symmetry is broken in the superconducting phase makes the 
physics of this model rich and complicated. Our model appears to be a holographic real- 
ization of the type of scenario described from a formal perturbative field theoretic point of 
view in ref. |T8] . Transport coefficients such as the speed of second sound and conductivities 
will depend on which direction we decide to look. Such a breaking of rotational invariance 
is not unheard of in real world materials. To pick a particularly simple example, a ferro- 
magnet will break rotational symmetry when the spins align. Ref. [9] emphasized a possible 
connection of this SU(2) model to a p-wave superconductor, where the order parameter for 
the phase transition is a vector. 

Of the real world materials that we considered, superfluid liquid helium-3 perhaps comes 
closest in approximating the physics of our model. Liquid helium-3 at very low temperatures 
is a p-wave superfluid. Two fermionic helium-3 atoms pair up to form a loosely bound 

^See refs. [T3[TH] for other nice analytic results for this class of models. 
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bosonic molecule with weak interaction between the orbital and spin degrees of freedom of 
the electrons [[I9j. The orbital and spin angular momenta are both equal to one, and the 
order parameter is often written Aai where a indexes the spin angular momentum and i 
the orbital angular momentum, in surprisingly close analogy with our j^. There are many 
stable phases of superfluid helium-3, depending on the pressure, temperature, and applied 
magnetic field. The A phases are known to break rotational symmetry. 

Despite plausible similarities between the symmetries of our model and various real 
world materials, there is one crucial difference. While the order parameters for these real 
world materials may have vector or tensor structure, they are not currents, and the signature 
of the phase transition is not the production of a persistent current. In contrast, our model 
has (jf) + 0. 

We begin in Section 2 with a discussion of the SU(2) model and the probe limit. We 
choose to work in a limit in which gravity is weak and the non-abelian field does not back 
react on the metric. Thus, at heart, in this paper we will be solving the classical SU(2) 
non-abelian Yang-Mills equations in a fixed background spacetime, that of a Schwarzschild 
black hole in jddS'sj^ 

In Section 3, we find a solution to the Yang-Mills equations near the phase transition. 
This power series solution in the order parameter and superfluid velocities allows us to 
demonstrate that the phase transition is second order and to calculate the speed of second 
sound from thermodynamic identities. 

In Section 4, we make some formal remarks about the current-current correlation func- 
tions for our model. We discuss the Ward identities that these Green's functions satisfy 
and some of their discrete symmetries. We also review how to calculate these two-point 
functions using the AdS/CFT correspondence. 

In Section 5, through a study of fluctuations about our solution near the phase transition, 
we extract the current-current correlation functions in the hydrodynamic limit. From the 
location of the poles, we independently confirm the speed of second sound calculated in 
Section 3. We are also able to calculate various damping coefficients and see explicitly that 
the Green's functions satisfy the non-abelian Ward identities. In the last part of Section 5, 
we consider the u; — > and ^ limits. From these limits we extract conductivities and 
also demonstrate that the system obeys a type of London equation. 

2 The Model 

Consider the following gravitational action for a non-abelian gauge field -F^^ with a cosmo- 
logical constant A: 

S=^l d'^+^x^g [R-2K)-^j d'^+^x^g F^^F'^^^ . (2.1) 
Our gauge field can be re-expressed in terms of a connection as follows: 

FIb = dAA% - dsA'X + f\,A\A% , (2.2) 

Attempts to solve the full set of coupled equations for a non-abelian black hole go back many years 
[501111]. See refs. [HllS] for reviews. 
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where /"^^ are the structure constants for our Lie algebra g with generators Tq such that 
[Ta,Tfe] = ifab^Tc- We win take g = su(2) where Ta = (Ta/2, aa are the Pauh spin matrices, 
and the structure constants are fabc = ^abc^ 

The equations of motion for the gauge field that follow from this action (2.1) are 
DaF'^^^ = which can be expanded as 

^^paAB ^ ja^^^h^pcAB ^ q ^ (2.3) 

Einstein's equations can be written 

Rab + (a - sab = ^ (2FSaF^^b - ^F^z.F'^^^^ab) • (2.4) 

A solution to these equations in the case of a negative cosmological constant, A = 
—d{d— l)/2, is a d+ 1-dimensional Reissner-Nordstrom black hole with anti-de Sitter space 
asymptoticsj^ The only non-zero component of the vector potential can be taken to b^ 

Al = + pu'^-'^ . (2.5) 

Thus we are using only a U(l) subgroup of the full SU(2) gauge symmetry; this black hole 
solution requires only an abelian gauge symmetry. The line element for this black hole 
solution has the form 



f{u)dt' + dx' + ^ ' 
where dx'^ = dx^ + dy^ + dz^ and the warp factor is 



ds^ = ^ 



(2.6) 



/ \ 2d-2 / \d 

fiu) = l + Q^-) -(1 + q2)M ^ (2.7) 

\UhJ \UhJ 



with the charge Q being defined as 



The horizon is located at u = u^, and the Hawking temperature is 

One subtle issue to be addressed is the boundary conditions for ylf at the horizon u = Uh 
and at the boundary n = of our asymptotically AdS space. At the horizon of the black 
hole, we must work in a local coordinate patch for the gauge potential such that has a 
well defined norm, < oo. Given the form of g**, we actually require that At{uh) = 0. 



The g = su(2) indices a,b,c, . . . are raised and lowered with the Kronecker delta 5^. The capital indices 
A, B,C . . . are raised and lowered with the five dimensional space time metric Qab- We will also shortly 
introduce Greek indices . . . which will be raised and lowered with the four dimensional Minkowski tensor 

77"" = (-+++)■ 

®We set the radius of curvature L = 1. 

'^The notation p is meant to evoke a charge density. The actual charge density p = — (d — 2)p/ . 
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Our gauge potential (2.5) is well defined globally, at both the horizon and the boundary, 
provided 

p=-^/nf-2. (2.10) 

The boundary value of the gauge potential, (0) = /i, is interpreted via the AdS/CFT 
dictionary as a chemical potential in the dual field theory. As such, fi is an external 
parameter of the field theory, and we should restrict our class of gauge transformations 
to those which do not affect /i. For example, the class of abelian gauge transformations 

— > + dtA where A = ct is ruled out by this restriction. 

Refs. jSl El [To] made the observation that below a critical temperature (or alternately 
above a critical chemical potential) , this charged black hole undergoes a second order phase 
transition. The component of the gauge field A^. develops a profile which not only spon- 
taneously breaks the remaining U{1) of the SU{2) symmetry but also breaks rotational 
invariance. These authors considered the case d = 4 and interpreted the dual 2+1 dimen- 
sional field theory as a p-wave thin film superconductor. More recently Basu et al. |13) 
looked at the d = 5 case where they interpreted the dual 3+1 dimensional field theory as a 
pion superfluidj^ 

In this paper, we shall make two simplifying assumptions. The first is to take the probe 
limit, as was done in [9l[ini[l3]. In this limit, K^/^f^ — > and the gauge field does not back 
react on the metric. The metric remains that of an uncharged black hole in anti-de Sitter 
space with warp factor 

f{n) = !-(-) . (2.11) 

Next, we restrict to the d = 5 case because of the observation made in ref. |tl3j that the 
zero mode inducing the phase transition has an analytic form. Given this form, we are able 
to compute analytically many properties of the field theory close to the phase transition 
including Green's functions, conductivities, diffusion constants, and the speed of second 
sound. We show explicitly that the phase transition is second order. 



3 Critical Behavior 

We specialize to AdS^ where we can construct an analytic solution to the gauge field 
equations of motion close to the phase transition. We give the solution as the first few 
terms in a power series in three small parameters: the order parameter e = g'^{j'f)/2 
and chemical-potential-like objects we call superfluid velocities, A^{u=0) = = v^^ and 
Ay{u=0) = Ay = v±. Velocity is a bit of a misnomer here as the objects and v±, like the 
chemical potential have mass dimension one. The name is motivated by their canonical 
conjugacy to the currents and j|. 

From the AdS/CFT dictionary, the currents (ja) and external field strength A'^ in the 
field theory can be determined from the small u expansion of the bulk gauge field 

^;i = -A^ + ^5'0>' + ... • (3.12) 

®See also [14] for a 3+1 dimensional system with similar symmetries and qualitative behavior but a more 
complicated action. 
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3.1 The Background 



We begin with small steps and construct the solution in the limit v± = = 0. In the probe 
approximation, the equations of motion for the gauge field take the form 

V,Af=^-^Al and V,AI = -^-^AI, (3.13) 

where we have defined the linear second order differential operators 



Vt = dl- -du and V, = Vy = dl+(^-^--]du. (3.14) 
u \ J u 



To keep the equations simple in what follows, we choose to put the horizon of the black 
hole at = 1. To restore units, dimensionful quantities such as the chemical potential //, 
frequencies lo, and wave- vectors k should be replaced with the dimensionless combinations 
ixuh-, ujuh, and kuh, respectively. 

As pointed out by ref. [1^, when A'l = 4(1 — u^) there is an analytic solution to the 



second equation of (3.13) that is regular at the horizon, of the form 



^ (1+^2^2 



(3.15) 



From eq. (3.12 ), the meaning of e in the dual field theory is, up to normalization, that of an 
expectation value for the non-abelian current {j^) = 2e/g'^. The existence of the solution 
(3.15) indicates that the superfluid phase transition occurs when /i = 4j^ Given this zero 



mode, we look for a general solution to eqs. (3.13) as a series expansion in e: 



^2 

= ey^^^ + eV + e V + O(e^) , (3.16) 
(1 + u^j^ 

The solution describes the system for > 4. Our strategy will be to fix the expectation 
value of (j^) = 2e/g'^ but to allow the chemical potential to be corrected order by order: 
/X = 4 + e'^dfii + e^5//2 + . ■ .. Thus, in solving the differential equations, we require the 
boundary condition that the 0{u'^) term in Wi vanish while (pi{0) is allowed to be nonzero. 
The differential equation governing (f)i is 

4n^ 

= , (3.18) 



which has the solution 



M 2^x , 1 /^r; 2 8n2(l + 3u2 + uS , 

(1 - u^)6,, + - [5u^ - \,^^.)3 ) ■ (3-19) 



^There are in fact a countable set of such zero modes with /i = 4fc where fc is a positive integer. We 
discuss these higher zero modes in Appendix [X] As the higher zero modes have higher free energy, they 
should not affect the phase diagram of the system. 
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We applied the boundary condition that 0i vanish at the horizon. Also, corresponds to 
a shift of the chemical potential by e^S^i. The value of is constrained by the solution 
for wi, as we now see. The differential equation for wi is 

„ l + Zu^ , 16 ^ ^ 

^1 71 4T^^l + n , 2\2 ^l = - J. 2\l-\ I 2^4 '^'^! " (3-20 

M(1 — M^j (1 — + M^j^ 

We require the boundary conditions that w\ be regular at the horizon and vanish at the 
boundary (u = 0). These conditions leave us with the solution 

cu^ u4(39m6 - SSlu"^ - 819n2 - 369) 13u2ln(l + u2) 
""^ ~ (1 + ^/2)2 + 20,160(1 + n2)5 ^ 1680(1 + u2)2 ' 

and the constraint 

= iio ■ 

The term in w\ proportional to c is just the zero mode, and, consistent with our strategy, 
we set c = 0. 

For the free energy calculation we perform below, we also need the next order corrections, 
and W2- The expressions are too cumbersome to repeat here. The structure and boundary 
conditions are analogous to the case of 0i and w\ considered above. 
The near boundary expansion of our solution takes the form 

A\ = tu^ + 0{u^) , (3.23) 

A = (4 + ^+5^26^ + 0(6^)) (3.24) 

/ 281e2 /1343 - 13651n2 , \ . 2 4n 

-(' + ^" 1 2,822,400 ^^2)e^ + 0(e^)j.^ + 0(u^), 

13(-4,015,679 + 5, 147,520 ln2) . ^ , 

= 75,866,112,000 " ^^"'^^ 

These expansions match well with numerical solutions that we found close to the transition 
temperature. 

3.2 Superfluid flow 

In this section, we generalize the background above to allow for the possibility of a superfluid 
flow. In terms of the bulk solution, this generalization requires turning on a constant value 
of Ay{u=0) = v± and A^{u=0) = W|| at the boundary corresponding to a non-zero superfluid 
velocity {v\i ,v±,0). The differential equations describing this background are a modification 



where 



of eqs. (3.13 ): 



u'fVxAi = g^^'^iAlAlAi-AlAlAl) , (3.26) 
u'fVxAl = g^^''{AlAlAl-AlAlAi) , (3.27) 
-^t^uAf — AfduAj- = fA^duA"^ — fA'^duA^ + fAyduAy — fAyduAy , (3.28) 

where we set Aj^ = A'^ = A'^ = 0. The repeated covariant A indices on the left hand side 
are not to be summed over. As before, we solve this system in a small e expansion, but 
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we also add another small expansion parameter ~ t;_L ~ f y . There is a non-uniformity in 
the limit v± ^ and 0, and we find two branches of solutions for small values of the 
superfluid velocity. In the case where v± > v^^, we find 

Aj = Oie'), 

(1+^2)2 ^ -L W 24(1 +^2)2 

- "'^;i(rT^ + ^ ii+^v 24(i+n2)2 +•••' 

A', = 4(l-n2) + i(.J + .i)(l-^2) 

2 + f J (1 _ n2)(71 + 3n2 - 627^^ - 279^6) 

6720(1 + ^2)3 + ■ ■ ■ ' 

3 _ _ ^vw vf\+^l M^(3 + 9n^+4n^) 
^ ~ """"^^^^^I 144(1 + m2)3 +•■■ ' 

3 _ _ n^(3 + 9^^ + 4^4) 

2^ ~ " ^ 144(1 + ^2)3 + ■ • • • 



In the case v± < f y , we find 



A\ = e . A . J +■■■ 



V\\ 4(1 +u2 

2 2^^*^(^^^ -21n(l +M^)) 

(1+^x2)2 + '^^^ + ^11 > 24(1+^.2)2 
U_L ,2 2n^-L '"^(^^ ~ 21n(l + u2)) 



" %|, (l + n2)2+ '^''^ + ^11^,1 24(1 + ^.2)2 + 



6 

+ (1 _ ^2)(7i ^ _ g27^4 _ 279^6) 
^^2 6720(1 + ^2)3 + • • • 



f^J + ^^i n2(3 + 9u2-2u 



At — v\\ e " , + • • • 



4\ 



2^3 



«ll 288(1 + u 
,3 ^ _ 2^± i + ^i^'(3 + 9u2-2u^) 

2/ ' ^;|| ^;|| 288(1 + U2)3 

These solutions can be used to compute the speed of second sound perpendicular and 
parallel to the order parameter A^.. In a two component fluid, there are typically two 
propagating collective modes, ordinary and second sound. In our probe approximation, we 
see only the superfluid component, and the single collective motion available to us we call 
second sound. From our holographic perspective, ordinary sound would involve fluctuations 
of the metric so it is suppressed in the limit /g'^ — > 0. 

The speed of second sound, like that of ordinary sound, can be computed from derivatives 
of the state variables. From ref. |llj . the second sound speed squared in this probe limit 



8 



should be 



dj /dv 
dp/djjL 



(3.41) 



D=0 



From eq. (3.12), the values of the charge current j = {jf) and the charge density p 



01) 



can be read off from the order pieces of Af and A'l, respectively. 

Because our system is not rotationally symmetric, the speed of second sound will depend 
on the direction of propagation. Let c± and cy be the speeds perpendicular and parallel to 
the order parameter A]., respectively. The speed c± can be computed from the background 
solution v± > I'll while can be computed from the solution with f|| > v±. In the case 
v± > v\\ = 0, we find that 



A' 



p + ^{840 - 281 fi)u'^ + ... , 

140 , 2 
v± - Y^vxip - 4)n + . . . . 



and hence, up to higher order corrections in e, 

^-^ ^ 13,488^ ^ 281 



(^-4) . 



(3.42) 
(3.43) 

(3.44) 



(We used the fact that ^ — 4 ~ 71e^/6720, which can be read off from (3.24).) In the case 
vu > f_L = 0, at leading order remains the same but now we need 



At 



70 



II - — U||(/i-4)n + 



1 71 



70 



We find that 

^11 ^ 2 13,488^ ^ 281 
We confirm these results for cm and c± in Sections 



(m-4) 



(3.45) 



(3.46) 



5.2 



and 5.3 through an analysis of the 



hydrodynamic poles in the current-current correlation functions. For numerical results valid 
when e is not necessarily small, see Figure [2j 

These perturbative solutions in v± and can also be used to analyze the phase diagram 
of the system near the critical point pc = 4. At the critical point, we expect the order 
parameter to vanish, so e = 0. The value of at u = can be reinterpreted as the value 
of the chemical potential. These two facts give us a relation between the chemical potential 
and superfluid velocity along the critical line separating the two phases. For superfiuid fiow 
parallel to the order parameter, we expect 



b " 

while for fiow perpendicular to the order parameter, we have instead 



(3.47) 



(3.48) 
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3.3 The Free Energy 

We compute the contribution to the free energy from the gauge field term in the on-shell 
action: 



_ _ aAB 

/3V0I3 du ( 1^ ^3^2 il \/Q \\\2 , 1 



^ {(d^A\f - m{d^A\f + (3.49) 



2 

For a background where = and 

= (4 + (5/xie2 + (5M2e^)(l - u^) , (3.50) 

the on-shell action is 

/3V0I3 , 71 2 51,145,217 4979 \ 4 



- V l^" ^ 1- 2.370.816.000 + i7(M(io ^ j + ' j ' P'"^'' 

Here V0I3 is the spatial volume of the field theory while (3 = 1 /T is the inverse temperature. 
For the background where A], 7^ has condensed, we find in contrast that 

S., = ^ (64 + + f - JWIL + J?!^ 2) + 0(.»)) . (3.52) 

4^2 \^ 210 V 2,370,816,000 176,400 J ^ ') ^ ' 

The difference in the values of the two on-shell actions is 

/3V0I3 / 71 , 



,3AP = Sv„ - S., = ^ + J . (3.53) 

Now AP can be interpreted also as a difference in the free energies because the free energy 
(in the grand canonical ensemble) is minus the value of the on-shell action. That AP < 
implies that the free energy of the superfluid phase is smaller and thus the superfluid is 
stable. 

Moreover, from the fact that the free energy difference scales as e^, we see that the phase 
transition is second order. For small e, ~ (/i — /ic)^- If we restore dimensions, then /x 
should be replaced by jiUh = fi/irT. Thus, ~ (Tc — r)^. The derivative of P with respect 
to temperature is continuous but non-differentiable. 



4 Formal Remarks about Green's Functions 

A field theory with a non-abelian global symmetry, such as SU(2), has by Noether's The- 
orem, a conserved current ja which transforms under the adjoint representation of this 
symmetry group. In this paper, we are interested in Green's functions for this current, in 
particular the Fourier transformed retarded current-current correlation functions: 

G'abiP) = i I d'xe-'P-mj^{x)jaO)])9{t) . (4.54) 
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If the symmetry is non- anomalous, then we can weakly gauge it by coupling the current to 
an external gauge field ^J^. Gauge invariance then implies that the correlation functions 
obey a series of Ward identities. For the one-point function, the covariant derivative of the 
current vanishes: 

0={d,dl + U^Al){ji!). (4.55) 

More usefully for the present discussion, there is also a Ward identity for the retarded 
two-point function. We give here the Fourier transformed version: 

= {ip.S'a + /afe^4)G^d (P) + f-'i'ifc) ■ (4.56) 

For our gravitational system beyond the phase transition, the gravitational bulk values 
of both and A], are non-zero. The AdS/CFT dictionary allows us to read A'^ and {ja) 



from the near boundary expansion of using the relation (3.12). For the system under 
consideration here, of the components of the external gauge field only A^ = fi is non-zero 
in the field theory. We have two non- vanishing components of the current, (jf) and (jg). 

Below, we compute the Green's functions that describe the response of the system to 
an external gauge field in the third isospin direction: G^^. The relevant Ward identities 
componentwise are 

= ip^G^r - (jT)'^a25'^% (4.57) 
= ip^G^3^ + /iGf3 + (jf)5-^ (4.58) 
= ip^Gt', - fiGt^ . (4.59) 

Our Green's functions below obey this set of Ward identities. 

Another important observation for the Green's functions under consideration is the 
symmetry under swapping the indices. We observe that 

G'atiP) = i-^)^'~^''^GZ{p) , (4.60) 

where (j){a,b) is equal to —1 if either a = 2 or 6 = 2, but not both, and 1 otherwise. This 
symmetry follows from the discrete symmetries of the system. Given that our currents 
are even under PT, i.e. parity and time reversal, if PT were a symmetry of the state, we 
would expect the Green's functions to be symmetric under an index swap. Our state is not 
symmetric under PT, but it is symmetric under PT times a Z2 operation on the 5u{2) Lie 
algebra, ai — > — ui and (T3 — CT3. 

4.1 Computation of two-point functions 



To compute the current-current correlators (4.54) in the probe approximation, we perturb 
the background gauge field by sending 

A'X^A'X + SA'X . (4.61) 

Consequently, the corresponding field strength F^^ changes to F^^ + f%^-, with /^^ given 

by 

fls = QaSAb - SbSAa + e'''''5A\A% + e''''^A\bA% . (4.62) 
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From (2.1 ), one can see that the quadratic action for SA'\ is 
which gives the hnearized equations of motion 



gfABf'''"' 



Jab + e ^ Jab + « oA 



^abc /lAb fc 



abcx A Ah r?c 



AB 



(4.63) 



(4.64) 



The quadratic action (4.631 is in fact not weh defined because the integrand diverges as Intt 



at small u as we will see. We will regulate this divergence using holographic renormalization 



For definiteness, we will only analyze the case where the background gauge field doesn't 
depend on t or x and where its radial components vanish. We choose a similar gauge 



for the perturbations by requiring 5A'^ = 0. Equations (4.64) can be solved approximately 
in the limit of small u. An appropriate series expansion in this limit is 



M^(t, X, u) = a'Ji.it, x) + a1{t, x)u^ In n + /3^(t, x)u^ + 



(4.65) 



for some vector- valued functions a(t, x), a(t,x), (3{t,x), etc. The values of a and /3 are 
the only ones that can be specified independently; all the other functions appearing in this 
expansion, namely a and higher order corrections, can be expressed in terms of a and (3. 



Plugging (4.65) into (4.64) and looking at the term with the lowest power of u in the 

s a relation between a'^ and a": 



equation with B = one finds a relation between a1 and a": 

1 



a,. 



(4.66) 



M = 



Upon integration by parts in (4.63), the unregularized on-shell quadratic action can be 
written as 

oon— shell 



25^ 



^ ^ (fx-{5A^^){du5Al) 



M=l/A 



(4.67) 



The divergence that arises as one takes A — > oo comes from the duSA'j^ term whose most 
divergent piece goes like u In u at small u. This divergence can be regulated by adding the 
counterterm 



ct 



In A 



d'^x 5A' 



u=l/A 



(4.68) 



Note that this counterterm depends only on the values of the gauge field on the surface u = 
1/A and on its derivatives along this surface, as required by holographic renormalization. 

As a side note, a simpler formula for a" can be found if only A^ and approach 
non-zero values at the boundary of AdS. In this case, only a'f, is non-zero and is given by 



1 



Assuming that a°(t,x) = a"e ^'^^ then 



1 



a,. 



if - u^')al 



(4.69) 



(4.70) 
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To compute the Fourier transformed two-point function, we first Fourier transform the 
regulated on-shell action 



^on-shell = 1 I ^ <{-p){P%p) + Calip)) , (4.71) 

where c is an arbitrary constant introduced by the regularization procedure. Although 
such an action is not a generating functional for the retarded Green's function, using the 
procedure outlined by Son and Starinets [25] , we can identify the retarded Green's function 
a£3 

2 d [(3;{p) + ca °(p)] 
dal{p) 

The linear response of a system to a perturbation a\,{p) is then a current density of the 
form 

{fM = • (4.73) 

For most physical questions, the ambiguity in the choice of c should be irrelevant. More 



precisely, one can see from (4.66) that schematically a = dda + da + a, so G"j^(p) is 
ambiguous up to an additive term analytic in p. Its Fourier transform G^{x) is ambiguous 
up to an additive term of the from ci(5^(x) + C2dx6^{x) + c^''d\dpS'^{x), where ci, C2, and 
Co'' are constants that depend on the particular Green's function we are computing. Since 



in position space equation (4.73) reads 



{j^{x)) = / d'x' GZix - x')al{x') , (4.74) 



it follows that the ambiguity in the choice of c does not affect the result of {ja{x)) if 
a^(x) = 0. In particular, the late-time, large-distance response of the system to localized 
sources is not affected by this ambiguity. There are many subtleties in these calculations. 



5 Fluctuations 

To calculate the Fourier transformed retarded current-current correlation functions, we need 
to study fluctuations of the SU(2) gauge fields A'^{x) in our black hole background. 

In the superfluid phase, the expectation value of the order parameter A], ^ breaks 
rotational symmetry and makes our task richer and more complicated than in the rotation- 
ally symmetric case where only ^ 0. In the rotationally symmetric case, it would be 
enough to consider a fluctuation with a time and space dependence of the form 6^*"^*+**^^. 
Given the breaking of rotational symmetry, we should in principle consider a more general 
dependence where we allow for motion both parallel and transverse to the order parameter: 
^-luit+ik^x+ikyy ^ Because of the complexity of the full result, we shall not present a full 
accounting of all the Green's functions here. Instead we will content ourselves by studying 
various informative limits where either kx = ^ 01 ky = {). 

^°For a more precise discussion of how to derive these Green's functions from an action principle and 
generating functional, see ref. |26) . See also ref. |27) and the discussion in Appendix C of ref. [28] . 
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We make a few other additional simpHfying restrictions. Following in the footsteps of 
refs. [9l [To] where the third isospin direction was interpreted as the U(l) of electricity and 
magnetism, we will consider Green's functions where at least one of the SU(2) isospin indices 
is equal to three. In other words, we are interested in the linear response of the system to 
external electric and magnetic fields. 

The last simplifying restriction is to limit our study to the hydrodynamic regime, where 
the order parameter, the frequency, and the wave-vector are small compared to the tem- 
perature. In our dimensionless notation, e,k,LO <^ 1. It is only in this limit that we have 
analytic results although it is straightforward to calculate the Green's functions numerically 
beyond this regime. 

We work out the Green's functions in five cases. The first and simplest case, for which 
we give the most detailed description of the calculation, is for a fluctuation transverse to 
the order parameter and a wave vector transverse to both the order parameter and the 
polarization of the fluctuation. We call this fluctuation the pure transverse mode. We next 
consider fluctuations that correspond to a second sound mode in two different limits, one 
where the sound is propagating parallel to the order parameter and one where the sound 
is propagating transverse. These two sets of fluctuations give us independent confirmation 



of the speeds of second sound computed in Section 3.2 from thermodynamics. Finally we 
consider fluctuations that correspond to a diffusive mode, again in two different limits, one 
where the diffusion is parallel to the order parameter, one in which the diffusion is transverse. 
(In Appendix pi we attempt to give a picture of the general case and how the different limits 



of these Green's functions fit together.) In Section 5.6 we discuss conductivities and the 
London equations. 

In what follows, to avoid cumbersome indices, we define new variables for the background 
values of the gauge field: 

Al = W and Af = ^. (5.75) 
5.1 Pure transverse mode 

The pure transverse mode is described by fluctuations of the field Ay with only z spatial 
dependence. We decompose the fluctuations into Fourier modes: 

SAliu, t, z) = ay{u)e-'^'+'^' . (5.76) 

These modes transverse to the order parameter A\ decouple from the other fluctuations of 
the gauge field and are governed by the differential equation: 



(fc2 + |^2)/_^2 

Vyay = ay , (5.77) 



where Vy was defined in eq. (3.14|. 



Near the horizon u = 1, we find that ay ~ (1 — n)^*"^/^ satisfies either ingoing or outgoing 
plane wave type boundary conditions. Consistent with the presence of an event horizon, it 
is natural to choose ingoing boundary conditions (the minus sign in the exponent). This 
choice leads to retarded, as opposed to advanced. Green's functions in the dual field theory 
At the boundary it = of AdS, we would like the freedom to set ay{Q) = ayo to some 
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arbitrary value of our choosing, corresponding to perturbing the dual field theory by a small 
external field strength. These two boundary conditions along with the differential equation 
uniquely specify the functional form of ay. 

While an analytic solution to eq. (5.77) does not appear to be available, one can easily 
solve this equation in the limit of small iv, k, and e. We can write the solution for ay, valid 
to order e^/c, e^w, /c^, and uj'^, in the form 



CLy — ClyO 



1 - 

1 +li2 



We find 



^yuje 



■^yk 



(l + e^aye + e^ujay^e + k'^CLyk + ^^'^CLy^j + ..) 



n2(3 + 9'u2 + 4ii'^) 

144(1 + ^2)3 ' 
iu2(12 + 27^2 + ISii^) 
864(1 + u2)3 



2 In (u) In 



l + n2 
l-n2 



(5.78) 

(5.79) 
(5.80) 
(5.81) 



The expression for ay^^ is too cumbersome to give here. Near the boundary, this solution 
(5.78) has the expansion 



lUO 



^ iioe^ a;^ln2 1 o i2\ 
ay = ayo + ayQ\ — - — - — — + -(w -k) 



\n{u) + ... (5.82) 



From this near boundary expansion and eq. (4.73), we can calculate the two-point function 
for the current in the hydrodynamic limit: 



2 f iu) iue^ In 2 



+ {u? - k^)c ]+.... 



(5.83) 



^2 2 48 72 4 

Note that the counter-term ambiguity, proportional to an arbitrary constant c, is of the 
form predicted in eq. (4.70). 



5.2 Transverse sound fluctuations 

In general, second sound modes are expected to produce poles in the density-density cor- 
relation function. We thus need to consider fluctuations in the conjugate field A^. If we 
consider sound modes moving transverse to the order parameter, we can take the fluctua- 
tions to have a y dependence but no x dependence. The self-consistent set of fluctuations 
to consider that couple to 6A'l(u,t,y) are 

dAl{u,t,y) = al{u)e-''^'+'^y , 

5Al{u,t,y) = al{u)e-'''^+^^y , (5.84) 

5Al{u,t,y) = aS(u)e— 

where a = 1,2. 
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The four fluctuations satisfy four second order ordinary differential equations and one 
first order constraint: 



Vyal 



Vtai 



4 + 



-UJ 



P 

P 

oj'^ + W^f o kuj o 

—a " 

y p 



P 



P 



2iuj^a]. 
ikW n 



/ 



at + y flx + 



/ 

2W^ 



f 



-a. 



lUJ 



1 



dual + ikdual + Wdual - [duW)al , 



(5.85 
(5.86 
(5.87 
(5.88 
(5.89 



where P^, Vx, and T>y were defined in eq. (3.14). We checked that the derivative of the 
constraint equation (5.89) with respect to u is a hnear combination of all five differential 
equations ( 5.85 )-( 5.89). Thus if a solution of the first four differential equations satisfies 
the constraint for some u, it will satisfy the constraint equation at all u. 

There are seven integration constants associated with this linear system ( 5.85 )-( 5.89). 
If we look at the horizon of the black hole at u = 1, we find seven different kinds of behavior. 
There exist six solutions that have plane wave behavior for a^, a^, and near the horizon 
of the form 



(1 



There is also a pure gauge solution, 



-lUO 



ik 



-w . 



(5.90) 



(5.91) 



As in the pure transverse case, we choose pure ingoing boundary conditions corresponding 
to (1 — n)"*'^/^ behavior. At the boundary n = of our asymptotically AdS space, we 
would like to be able to perturb the system with arbitrary boundary values of af and 
but set the "unphysical" components of the gauge field and to zero. These are four 
constraints and we have only three ingoing solutions. Thus we will also need to make use 
of the pure gauge solution to enforce our n = boundary conditions. 

We solved the system perturbatively in iv, k, and e. We present the results here in the 



limit where uj ~ 
form 



k' 



e . 



The near boundary expansion (u = 0) of the solution takes the 



{atpk + aypuj) 
V 



70ke (48A;2 + Se^ - 248iuj) 



atok + ayou^ ^u:e ^ _ ^^^^^^^^^ ^2 ^ ^^^^^2 ^ 



V k 

{atok + ayouj) 



V 



tj( 1120A;^ + 3A;2(117e2 - 1120icj) 



+ ^{e^ - 24ia;)(843e2 - 72,800ia;)) 
48 



[atok + ayoui) 



V 



k[ 1120/c'' + 3A;2(117e2 - 1120ia;) 



(5.92) 
(5.93) 
(5.94) 

(5.95) 



+ ^{e^ - 24ia;)(843e2 - 72,800ia;))u^ + 
48 / 
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Note, the expression {atok+ayQUj) is not homogeneous in our scahng hmit. We have included 
the leading corrections proportional to ato and UyQ. There are terms in the expansion 
proportional to u^lnu but they are subleading in uj, k, and e. 
The pole in this limit takes the form 

7i 

V = -72,800w^ + (43,120A;2 + 843e^)a;2 + — (4800A;'^ + 553A;2e2)w + 
V 6 

71 k'^f^ 

_ ^^^^4^2 _ ^^20A;6 + . . . . (5.96) 

16 

Let us study this cubic polynomial in lo in two different limits. First, if /c <C e, we find three 
poles with the asymptotic form 

, / 71 , 147,217iA;2 
" = 13:488^'- 947,532 +■•• ' ^'-'^^ 

843ie2 4,335,443i/c2 
~ ~ 72,800 ~ 15,397,395 + ••■ ' • 

The first two poles are propagating modes that we identify with second sound. Indeed, the 



speed of second sound agrees with the earlier result (3.44) from Section 3.2 The position 
of the third pole in this limit is determined mostly by the size of the order parameter e and 
so we associate it with the zero mode that causes the phase transition from the superfluid 
phase back to the normal phase. 

In the opposite limit, k ^ e, where the order parameter is small, the behavior should 
be close to that of the normal fluid. In this limit, we find 

/±ll-3i\,2 /±260,803 - 131,519A o 
= 26.644.800 f^-"'" 

ik"^ bie^ 

UJ = + . . . . (5.100 

2 2928 ^ ' 

The first two poles are associated with the zero modes that cause the phase transition from 
the normal phase to the superfluid phase and were discussed in ref. [9j while the third pole 
is associated to the diffusive mode of our conserved charge density. Indeed, the location of 
this diffusive pole is determined by the dynamics of the normal phase and was calculated, 
without the order correction, long ago in ref. [29]. As we vary e and k the number of 
poles cannot change. The two zero mode poles evolve into the sound poles of the previous 
limit while the diffusive pole becomes the zero mode pole of the previous limit. 

From these small u expansions, we can read off the eight Green's functions G'J^g, G^g, 



Gfg, Ggf, G3I, G^g, and G33. From the discrete symmetries (4.60), we can also read 



off four more Green's functions with the indices swapped. Note the prefactor atok + ayouj 



in the small u expansion. This structure is necessary to satisfy the Ward identities (4.57). 



As a further check, we consider a particular static limit of the density-density correlation 



function. From eqs. (4.73) and (5.95), we can read off the Green's function, 

2 k^ ( 1 \ 

G% = — ^— ( 1120/c^ + 3A;2(117e2 - 1120iw) + — (e^ - 2^iiS){M'i(? - 72,800za;) j . 

(5.101) 
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We are interested in the long wave-length limit of this Green's function: 

This long wave-length limit is equal to a thermodynamic susceptibility, 

lin,GS(O..) = = |. 



Given this relation, we see that eq. (3.42) agrees with eq. (5.102). 



(5.102) 



(5.103) 



5.3 Longitudinal sound fluctuations 

Longitudinal sound modes correspond to the case where the fluctuations in depend only 
on X. A self-consistent set of perturbations in this case is given by 

6A^{u,t,x) = a?(u)e-^^*+^'^^ , 

6Al{u,t,x) = a^(n)e-*^*+^^% (5.104) 

where a,b = 1, 2, 3. These fields satisfy the following six second order equations and three 
constraints: 

1 



Vtal 



f 
1 



- {2ikWal + (A;2 + W"^) + iWujal + kual + ik^al) , 
((fc2 + W^) 4 - 2ikWal + kujal - iWcoal + 2W<^al) 



f 
1 

P 
1 

P 
1 

P 

—^'af + iujdua] + ^duOt + ifkdual. , 

fW'al + <^'al - <^dual + icodua'^t + ifkOuul - fWduul , 

-fW'al + iujdual + fWdual + ifkdual , 



-2W<^al + ik^af - kuja\ + 2i^ujal - (^>^ + uj"^) al) 
-iWujal - kujaf - ik<^a\ - (^>^ + w^) - 2i<I>a;a^) 
-kual + iWuja} + W^a] - u?al) , 



(5.105) 
(5.106) 
(5.107) 
(5.108) 
(5.109) 

(5.110) 

(5.111) 
(5.112) 
(5.113) 



with T>t and T>x as defined in (3.14). Again, the three constraint equations are consistent 



with the second order equations in the sense that if they hold at some u, they hold at all u. 



The system (5.105 )-(5. 113 ) has nine integration constants. The nine possible behaviors 



at the horizon are of two types: six plane wave solutions for a^, a^, and that behave as 

(1 - u)^''^'^ (5.114) 
close to M = 1, and three pure gauge solutions given by 



ika 



al = ika^ - Wa^ 



a% = Wa + ika 



3 



(5.115) 
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where a"" are arbitrary constants. As in the previous sections, we require no outgoing modes 
at the horizon, which amounts to specifying three of the nine integration constants. The 
other six integration constants are specified in terms of the values of the fields at u = 0. 
In order to examine fluctuations in af and a^, we set their boundary values to Ofo and Uxo-, 
respectively, and the boundary values of the other four fields to zero. 

Solving the system (5.105 )-(5. 113) perturbatively in uj, k, and e under the scaling as- 
sumption a; ~ /c^ ~ e'^ , we find that the boundary behavior of the fluctuations is 



V 4 



(20,160fc^ + 843e^ - 72,800iw)u^ + . . . , (5.116) 



2 Ofofc + aa;ow 35A:2e 2 , 0-2 , oon ^ 2 , 

= — ^ — [Amk + 3ze -|-320a;)n + . . . , (5.117) 

al = + ""°^ 35A:6(48fc^ + 36^ - 320i^)n^ + . . . , (5.118) 

al = — ^ —(20,160fc2 + 843e2- 72,800iw)'u2 + a^o^w +--- ,(5.119 

V k ^ k 

a\ = ato + t,""""^ - (26,880fc" + 8436^ + 192fc^(79e^ - 840iu;) + (5.120) 

V 96 

-113,264ie2u; - 3,494,400u;2)n2 + . . . , 

al = a^o - t,""""" - (26,880fe^ + 843e^ + 192fc^(79e^ - 840za;) (5.121) 

V 96 

-113,264^^0; - 3,494,400w2)n2 + . . . . 
The pole here is again a cubic polynomial in uj: 

V = -72,800^^3 + (39,760^^ + 843e'^W + ^^(2688/0"^ + 451k'^e'^)uj 

6 

-— fc^e^ _ 53A;4g2 _ 280/^6 . (5.122) 
32 

We consider the roots of the polynomial first in the limit k <^ e: 



I 71 , 103,535 ,2 . X 

u = ±J ek ik^ + ... , (5.123) 

V 2 13,488 947,532 ' ^ ' 

843 .2 5,044,459 .,2 , 

a; = ^e + ... . (5.124 

72,800 15,397,395 ^ ' 

The first pair of poles correspond to second sound propagating in the direction parallel to 
the order parameter with a speed consistent with our earlier result (3.46). The third pole is 
related to the zero mode that causes a phase transition from the superfluid to the normal 
phase. Next we consider the limit k ^ e: 

±ll-3z ,^2 , ±192,553 - 95,119i 2 

UJ = k H e , 5.125 

130 26,644,800 ' ^ ' 

Lo = h . . . . (5.126) 

2 2928 ^ ' 

The two sound poles have evolved into the zero mode poles, while the zero mode pole has 
evolved into a diffusive pole. 
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From the small u expansion, we can read off a large number of Green's functions which 
we shall not bother to list. Similar to the transverse sound case considered above, the 



prefactor {atok + axo<^) in the expansion means that the Ward identities (4.57) will be 
satisfied. However, there is more structure here. Note that ika^ = Aa} — axo^v? and 
ika\ = — 4aj. In our hydrodynamic limit at leading order in uj and fc, these two equations 



are the Ward identities (4.58) and (4.59), respectively. 



Before moving on, we note that 

2 281 



lim^G**3(0,fc) = ^^ , (5.127) 



which agrees with eq. (5.102), but k here is parallel rather than transverse to the order 
parameter. 

5.4 Transverse diffusive mode 

In addition to the sound mode found above, in the limit A: <C e, we expect to find a diffusive 
mode in the current-current correlator. We begin with the slightly simpler case of a mode 
polarized transverse to the order parameter but propagating parallel to it, and follow in the 
next section with a mode polarized longitudinal to the order parameter but propagating 
transversely. Thus first we look for fluctuations in 5Ay{u^ t, x) and any other modes that 
couple to it. A self-consistent set of fluctuations to consider is 

M^(n, t, x) = a;;(u)e-*^*+^'=^ , (5.128) 

where a = 1, 2, 3. 

This set of fluctuating modes gives rise to the three differential equations at linear order: 



-D^al = ^ al-—al, (5.129) 



(P + iy2)^_^2_^2 2iW ^ 2iu:^ ^ 
T^xay = ^2 ay + —ay - -j^ay , (5.130) 

T>xay = j2 % + -J^% ■ (5.131) 

As before, we solve this set of equations perturbatively in the limit a; ~ /c^ ~ e'^. The 
small u expansion of the solutions, from which we may read off the Green's functions, takes 
the form: 

al = "^22keu;u'^ + . . . , (5.132) 
al = '^2ke{2ik^ + 2,uj)u^ + . . . , (5.133) 

al = ayo-"^-^(-1680k^ + 12k'^{29e'^ + 700iuj) -Gk'^ie'^ + 28ie'^LO - 10,780uj'^) + 
^ V 3360 V 

-hu;(9ie^-h 47666^0; - 109,200ia;2)^n2-F^ayoA;^M^ In n-h... . (5.134) 

The poles at leading order in this perturbative expansion come from a quadratic polynomial 
in ijj: 

2 ^ 3z(140fc2 + 3e2) (70P + 3e^)fe^ 
P = 65u; ^ w ^ . (5.135) 
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As before, we consider the roots of this polynomial in two limits. First we consider k <ti e, 
in which case we find 

2ik^ 

uj = ^ + -- - • (5-137) 

The first pole is associated with the zero mode that causes the phase transition from the 
superfluid phase to the normal phase while the second pole comes from a diffusive mode of 
the charge density. 

Next, we consider the limit k ^ e where we recover the zero modes of the normal phase, 

±11 -3i,2 ±33 -9i 2 

co = k^ -\ e^±.... (5.138) 

65 9100 ^ ' 



At leading order, the location of the pole is the same as that of eq. (5.99). However, the 
subleading order corrections are different. 

5.5 Longitudinal diffusive mode 

We continue the discussion by looking at modes polarized longitudinal to the order pa- 
rameter but propagating transversely. We consider fluctuations 6A^{u,t,y) and all others 
coupled to it: 

6Al{u,t,y) = a3.(n)e— *+^^^ 

6A1{u,t,y) = a?(n)e-*-*+*^^ (5.139) 
5Al{u,t,y) = a^(n)e-*-*+*'=^ 

where a,b = 1,2. This set of fluctuations obeys the five second order equations and two 
first order constraints: 



Vtal = 


y [-W^al - ik^al + k'^a] + kujal) , 


(5.140) 


Vta] = 


J (iWujal + {k^ + W^) at + kuja^ + ik^al) , 


(5.141) 


V^al = 


{ik^af + 2i^ujal - kojat - (^>^ ± u?) ay) , 


(5.142) 


V^al = 


{-ifkWal - kujal + {fW^ - $2 _ ^2^ ^2 _ -^^^^ _ 2i<!>ujay) 


,(5.143) 


T^xal = 


^ {{fk^ - Lo^) al + iWuja} + ifkWal + W'^aj) , 


(5.144) 


= 


-^'af + iujdua] + ^duaf + i/kdua^ , 


(5.145) 



= fW'al + <^'al - <^dual + iujduaj - fWdual + ifkdual . (5.146) 

At the horizon, there are two pure gauge solutions, three ingoing solutions, and three 
outgoing solutions. We discard the outgoing solutions and use the remaining degrees of 
freedom to choose the boundary values of the five fluctuations. In particular, we set the 
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boundary values of all the fluctuations to zero save for a^, which we set to a^o- The near 
boundary expansion of the solution takes the form: 

a] = ^^(fc2-3z..)u2 + a.o|n2 + ... , (5.147) 

at = -7p ^ u +... , (5.148) 

a\ = -"^sskeuj u'^ + . . . , (5.149) 

al = -"^Sikeik'^ - 3iuj)u'^ + . . . , (5.150) 



= a^^ + ^-^(840k^ -16k'^(13e'^ + 420ioj) -uj(e^ -48iuj)(9ie'^ +45b0oj) 
V 3360 V 

+?,k^{e^ + I2bie^uj-2,%,imuj'^)^u^ + ]^aMk^u^\nu + ... . (5.151) 
The pole, similar to the case considered previously, is a quadratic polynomial in oj: 

,2 , / aAu2 , 9«e^\ , , 3 22 



V = 130u;^ + 6ir + a; k^e^ - k^ . (5.152) 

35/35 



In the limit k <^ e we find a zero mode and a diffusive mode: 



u = -^e^ + ^k'' + ... , (5.153) 

4550 195 ' ^ ^ 

a; = -^ + ... . (5.154) 

In the opposite limit, we find two zero modes: 

±11 -3i,2 ±33 -9i 2 

Lo = e^±.... (5.155) 

130 9100 ^ 

The structure of the small u expansion of the gauge fields is again related to the Ward 

identities. We see that ika^ = Aaj — axo^v? and ikal = — 4a^, which are restatements of 



the Ward identities (4.58) and (4.59), respectively. 

We have thus far considered the current-current correlation functions in five special 
cases. While the results are simpler, our presentation has the disadvantage of obscuring the 
relationship between the various cases. We make some remarks about and provide some 
results for the general case in Appendix [B| 

5.6 Conductivity and London Equations 

In this section, we begin by studying the response of the system to a homogeneous, time 
dependent electric field, 5^ ~ e~"^^, and end with a discussion of the London equations. 
A homogeneous electric field should produce a current in the system via Ohm's Law. To 
investigate the conductivity in this long wavelength limit, we set k = for the two-point 
functions computed above. 

The case of an electric field orthogonal to the order parameter is simple; a current and 



nothing more is produced. From the pure transverse mode and eq. (5.83), we have 
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Reassuringly, this result agrees with the A; — > limit of the Green's functions associated to 
transverse sound propagation and the transverse diffusive mode. 

For an electric field parallel to the order parameter, the physics is richer. We find a 
current in the x direction but also oscillating (or precessing) charge densities associated 
with the one and two isospin directions: 

al = + 0x0 (^-^ + « Q + + • • • , (5.157) 

a] = axo^u^ + ... , (5.158) 

2 '^^^ 2 

at = -axo^u + ... . (5.159) 
lb 

This near boundary expansion agrees with the A; ^ limit of the expansions for longitudinal 
sound and diffusion considered above. The associated Green's functions are 

g2 \^ 96 \2 288 J 



G^3» = {-- + i{- + — ]o,]+... , (5.160) 



Identifying the electric field Ej = icuSAj and recalling Ohm's Law, the conductivities 



are related via eq. (4.73) to the retarded Green's functions. 



aUu;) = ^^ and ayy{u;) = . (5.163) 

The terms proportional to in Gff and G33 thus produce a pole in the imaginary part of 
the respective conductivities. As discussed in refs. |2i |30j, by the Kramers-Kronig relations 
(or by properly regularizing the pole) there must be a delta function in the real part of 
the conductivity, indicating the material loses all resistance to DC currents and suggesting 
the phase transition is to a superconducting state. While in refs. [21 130], the pole was seen 
only numerically, here we can calculate the strength of the pole analytically close the phase 
transition. Its residue is given by 



ReSaj=o(Txx = ^7j^ + • • ■ Res^=oo-ro = + ■ - - - (5.164) 



2 ie"^ _ 2 ie^ 

^— + ... iies^=0Tyy- 

In Figure [T] we show a comparison between numerical computations of the residues of the 
poles at w = in axx and ayy, along with the analytic approximation (5.164) close to 
T = Tc. 

In the Drude model for an ideal metal, the conductivity takes the form a = ip/muo where 
p is the charge density and m is the mass of the charge carrier. In the superconductivity 
literature (see for example |3lJ ) , the pole in the imaginary part of the conductivity is thus 
often related to a superfluid density. Because our system is not rotationally symmetric, the 
density to mass ratio defined in this way will depend on the orientation of the superfluid 
velocity with respect to the order parameter. The proper way to interpret this situation is 
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probably that a suitably defined effective mass of the superfluid depends on the direction 
of propagation. 

An important observation is that in our system, the u; — > and A; ^ limits of the 
Green's functions commute. The residue of the pole in the conductivity is related to the 
long wavelength limit of the current-current correlation function in the following way: 



iReSa}=o(Tjj = lim limGii{u,k) . 
The limit in the opposite order is related to a thermodynamic susceptibility: 

2" ' 



(5.165) 



and lim GgKO, h 



dv 



(5.166) 



where and v±_ are superfluid velocities It follows from eqs. (3.43) and (3.45) that 
d'^P djii 2 e2 , d'^P dj± 



dv\\ 



1^ 
'^2 96 



while 



dv± 



1^ 
'g^ 48 



(5.167) 



When combined with eq. (5.163), these results confirm eq. (5.164). To see how these limits 



commute in greater detail, the reader is referred to Appendix |B| and eqs. (B.179 )-(B.184 ) 
that give the general Green's functions in the inviscid limit kx,ky <^ e. 

As emphasized in this context in ref. [30], that the limits commute implies the system 
really does become a superconductor below T^.. Given that the limits commute, the system 
obeys a London type equation for small k and to: 



1^1 
'5296 



and 



52 48 y 



(5.168) 



If we now imagine the U(l) subgroup generated by T"^ G su(2) is weakly gauged, these 
London equations imply not only infinite DC conductivity but also a Meissner effect with 
London penetration depths that scale as A_l ~ An 1/e ~ (Tc — T)~^/2_ 



6 Discussion 



One of the nicest features of our results is their analytic nature. We were able to confirm 
a number of previous numeric observations [ll [HI [9l [TTl [32l [33] of this superfluid phase 
transition. In particular, we saw explicitly that the phase transition was second order; the 
difference in free energy between the phases scaled as {Tc — T)2 below the phase transition. 
We saw the order parameter grew as e ~ {T^ — T)^/2 below Tc and thus has a mean field 
critical exponent. We calculated the speed of second sound near the phase transition and 
observed that it vanished linearly with the reduced temperature c± ~ cy ~ {Tc — T). We 
looked at the pole at w = in the imaginary part of the conductivity and saw the same 

^^Note that to produce a perturbing magnetic field, we require a k that is transverse to the polarization 
of the current-current correlation function. A perturbation of the form 5A^ ~ e'*^ or SAy ~ e*''^ is gauge 
equivalent to zero and does not produce a response from the system. The Green's function in this limit 
vanishes. 
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Figure 1: Plots of numerical results for ^Res, 



olmaxx and ^KeSuj=olTaaxx as functions 
of temperature (solid lines), as well as analytical approximations at small e (dotted lines) 
given by equations (5.164). 



scaling, a ~ i{Tc — T)/lu, that had been observed numerically in a related model [2] and 
confirmed that the London penetration depth scales as A ~ l/(Tc — T)^/^. This laundry list 
of scalings (or critical exponents) is the same observed in the mean-field Landau-Ginzburg 
model of a superconductor. 

Close to Tc, eqs. (3.44) and (3.46) indicate that = 2cy, so one may wonder whether 
such a formula is valid away from Tc as well. Numerical evaluations show that this is not the 
case: see Figure |2j At small temperatures, approaches 1/3. Our numerical evaluations 
are not sufficiently reliable at small temperatures to see whether cf, has the same limit 
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There are some other results in this paper that are worth emphasizing. For superfiuTd 
velocities that are not too large, we were able to determine analytically the critical line in 
the temperature-superfiuid velocity plane separating the normal phase from the superfiuid 
phase. We calculated a large number of current-current correlation functions in the hy- 
drodynamic limit and verified that they satisfied the non-abelian Ward identities. We also 
investigated how the hydrodynamic poles in these correlation functions move around as a 
function of k and e. 

Optimistically, we hope that someday this system will be more than a toy model. The 
introduction described possible similarities of this system to helium-3 and p-wave super- 
conductors. Here we add a speculation about a possible connection with QCD. The SU(2) 
global symmetry of our model could be thought of as the residual approximate isospin sym- 
metry of QCD at low energies and our chemical potential an isospin chemical potential. 
The phase structure of QCD at non-zero isospin chemical potential has been discussed by 



^^In the case of a scalar order parameter and a phase transition that does not break rotational symmetry, 
we expect the speed of second sound to approach (d — as T — » 0. This limit follows from eq. (3.41 1 

and two observations: 1) At T = the Lorentz symmetry breaking due to the temperature disappears and 
the pressure can depend on /i and v only as P{i-i^ — v'^). 2) By dimensional analysis, when T = v = 0, 
P ~ /i'*. We would like to thank Amos Yarom for discussion on this point. 
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Figure 2: The squared speeds of second sound Cm and c\ as functions of the reduced tem- 



perature T/Tc (sohd Unes) as weh as analytical approximations given by eqs. (3.44) and 



(3.46) close to Tc (dotted lines). 



ref. [33]. Alas, there is no persistent current in the stable phases they discuss. 

In the future, there are some questions that we would like to address. While we focused 
on current-current correlation functions in the third isospin direction in this paper, it would 
be good to study the full set of Green's functions more carefully. Because of the charac- 
teristic magnetic properties of superconductors, it would be interesting to investigate the 
dependence of the correlation functions on an external magnetic field. 

Another interesting direction to pursue is the connection between this work and the 
membrane paradigm [35] where the horizon of a black hole, rather than the boundary of an 
asymptotically anti-de Sitter space, is thought of as a fluid. Related to this direction is the 
observation of ref. [8] that the fraction of the total charge density outside the black hole 
horizon scales as Tc — T close to Tc, suggesting that this quantity might be related to the 
superfiuid density. We would like to know if this analogy can be made more precise. 
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A Higher zero modes 



One can also find other analytical solutions to (3.13) when the chemical potential is = 
Ak, where k is an integer, corresponding to higher zero modes of A\. In fact, numerical 
exploration shows that these are the only values of /u for which zero modes occur. The 
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corresponding backgrounds are probably unstable, because they have higher free energy. 
When fi = 4k, the corresponding zero mode can be written as 



where e is an arbitrary constant and Pk{v?) is a polynomial of degree 2{k — 1) in v?'. It is 
not hard to find a recursion formula for the coefficients of Pk{x). Writing 

Pfc(x) = ^aWx- (A.170) 



and using equation (3.13), it follows that 

n{n + l)a(f ) + Ak{k - n)a^^}_^ - [Ak^ - 2k{2n - 1) + n{n - 1)] ajf^a = • (A.171) 

(k) 

The boundary conditions needed to solve this recursion relation can be taken to be a_( = 

(k) 

and = 1, the latter representing the normalization of the series. One can check that for 
each integer value of k the series terminates at order n = 2k, the last non-zero coefficient 



being a2(fc-i)- tablejljwe provide analytic expressions for the first few Pfc( 



Table 1: Formulas for Pk{u'^) with k < 5. 



k 


Pki 






1 


1 






2 


1 - 






3 


1 - 


12^x2 + f u4 - 


- 12^6 + 


4 


1 - 


24u^ + 135u^ 


- 240^6 + 135n8 - 24u^^ + 


5 


1 - 


40^2 + 412n^ 





B The General Case 

Consider a general fluctuation of the form 6A'^ = ^"(^i^e"*'^*"''*^^^"'"*'^^'^ for a = 1,2,3 and 
fi = t,x,y, in the limit where to ~ k^e ~ kyC. Solving the system of differential equations 
in the hydrodynamic limit leads to a set of Green's functions with a pole that is a flfth 
order polynomial in uj. Three of the roots of this polynomial correspond to the roots of the 
third order polynomial we found for the sound modes. Two of the roots of this polynomial 
correspond to the roots of the second order polynomial we found for the diffusive modes. 

In other words, the six modes we found in the longitudinal and transverse sound cases 
are actually different limits {ky ^ and k^ ^ 0, respectively) of three of the roots of 
this flfth order polynomial in iv. Similarly, the four modes we found in the longitudinal 
and transverse diffusion cases are different limits of two of the roots of this flfth order 
polynomial. 
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In the limit kx,ky <^ e, we find that the roots are 



CO 



LO 



UJ 



71 /l7^ ~ i(103,535A:2 + 147,217 
13,488 V 2 ^ ^.^^^^ 



-{2kl + kl) + 



947,532 



ne 



2 5g 

+ —i{2kl + kl) + 



4550 195 



843ie2 7i(720,637fe2 + 619,349 



72,800 15,397,395 
In the limit kx,ky ^ e, we find in contrast that 



+ ... . 
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(B.175) 

(B.176) 



+ 



33 



"9100 9100 



ki+{± 



260,803 



131,519? 



"26,644,800 26,644,800 



+ 



±— - —\ [kl + kl) + 

130 i3oy ^ y' 

2 ^ ^ 192,553 



(B.177) 



^•2 _|_ Z;-2 



95,119f 



LO 



"26,644,800 26,644,800 

U2 



33 



9i \ 



ki + [± 
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i(13A;^ + 5A;2) 
2928(^2 + k^) 
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(B.178) 



The case kx,ky <C e admits an inviscid limit where the two-point functions of the currents 
in the third isospin direction become 
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(B.179) 
(B.180) 
(B.181) 
(B.182) 
(B.183) 
(B.184) 



From eqs. (B.179 )-(B. 184 ) one can easily see that the limits w ^ and A; ^ commute. 



implying that there is a Meissner effect; see Section 5.6 
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Another interesting limit to take of this fifth order polynomial in uj is the limit w — > 0: 



7{) kl + kl j y ' y ' 1120 kl + kl J • 
From this limit, one can extract correlation lengths that scale as 1/e ~ (Tc — T) 



-1/2 
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